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Finite  Difference 
Interface 

Ordinary  Difference  Equation 


Reference  Wavenumber 
Spl i t-Step 
Wide  Angle 
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This  document  contains  the  slide  presentation  entitled  "The  State-of- 
the-Art  Parabolic  Equation  Approximation  as  Applied  to  Underwater  Acoustic 
Propagation  With  Discussions  on  Intensive  Computations,"  given  at  the 
108th  meeting  of  the  Acoustical  Society  of  America,  8-12  October  1984,  in 
.Minneapolis,  Minnesota. 


^  The  Parabolic  Equation  (PE)  has  applications  in  many  different 
scientific  fields  such  as  electromagnetics,  optics  theory,  quantum 


t.t  «, 


20.  Continued; 


'mechanics,  plasma  physics,  seismology,  underwater  acoustics,  and  others, 

The  subject  of  this  presentation  rs  PE  approximation  as  applied  to  underwater 
acoustic  wave  propagation.  A  review  will  be  given  on  past  contributions,  re* 
cent  developments  will  be  highlighted,  and,  looking  ahead,  we  will  discuss 
what  the  PE  method  can  do  in  order  to  stimulate  future  research  and  develop¬ 
ment,  as  well  as  applications.  Intensive  computations  with  respect  to  the  PE 
implementation  will  also  be  discussed, 
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THE  STATE-OF-THE-ART  PARABOLIC  EQUATION  APPROXIMATION 
AS  APPLIED  TO  UNDERWATER  ACOUSTIC  PROPAGATION  WITH 
OISCUSSIONS  ON  INTENSIVE  COMPUTATIONS 

The  Parabolic  Equation  (PE)  approximation  was  first  Introduced  by 
12  3 

Tappert  *  *  over  a  decade  ago.  Tappert' s  three  papers,  which  are 
outlined  In  vugraph  1,  are  often  referenced. 


Computer  Simulation  of  Long-Range  Ocean  Acoustic 
Propagation  Using  the  Parabolic  Equation  Method 

F.  0.  Tappert  and  R.  H.  Hardin 

Eighth  International  Congress  on  Acoustics, 
p.  452,  London  (1974) 


The  Parabolic  Equation  Approximation  Method  in 
■  Wave  Propagation  A  Underwater  Acoustics 

F.  D.  Tappert 

ed.  J.  B.  Keller  and  J.  S.  Papadakis 

Lecture  Notes  in  Physics,  Sprlnger-Verlag  (1977) 

Applications  of  the  Split-Step  Fourier  Method  to  the  Numerical 
Solution  of  Nonlinear  and  Variable  Coefficient  Wave  Equation 

R.  H.  Hard|n  and  F.  D.  Tappert 

SIAM  Review  15,  p.243  (1973) 

VUGRAPH  1 

In  recent  years,  many  Improvements  were  made  on  the  PE  technique  with 
respect  to  approximation.  Implementation,  and  application.  Let  us  begin  this 
discussion  by  a  brief  review  of  how  the  parabolic  wave  equation  was  derived 
by  Tappert. 

Consider  the  two-dimensional  Helmholtz  equation  In  cylindrical 
coordinates,  that  Is, 


i*  ♦  1 1*  ♦  »!i  ♦ 
17  r»r  7? 


kp  n2(r,z)  4-0. 


where  a(r, z)  Is  the  wave  field,  kQ  Is  the  reference  wavenumber,  n(r,z) 

Is  the  Index  of  refraction,  r  Indicates  the  range  direction,  and  z  Indicates 
the  depth  direction. 


les 
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k2(n2(r,z)  -  l)  u  -  0. 


(5) 
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Considering  only  the  outgoing  wave  In  the  range  direction,  we  see  that 
the  solution  of  Eq.  (4)  Is  the  zeroth  order  Hankel  function  of  the  first 
kind,  M^(k.r), 

Applying  the  farfleld  approximation,  kQr»l ,  to  the  argument  of 
H^V).  we  find  that 


(6) 


Using  Eq.  (6)  to  simplify  the  coefficient  (1/r  +  (2/v)  vp),  In  Eq.  5., 
we  find 


urr  *  (21ko>  ur  *  u«  +  ko("2(r’2)  -  *)  «  ■  0  • 

Dropping  urr,  based  on  the  paraxial  approximation, 
l«rrl«|21kourl,  we  find 


7  kQ(n2(r,z)  -  l)  u 
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(7) 


(8) 


This  Is  the  first  parabolic  wave  equation  derived  by  Tappert  and  Is 
referred  to  as  the  STANOARO  PE. 

I  would  like  to  take  another  approach  to  Eq.  (8),  which  you  will  find 
useful  In  observing  some  Important  physical  properties. 


'» ’**»-*»  -V '■ 
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(s '  \  -  % 


* ("2 " 11  +7IIyu  *  °' (9) 


Mote  that  Eqs.  (7)  and  (9)  are  the  same  If  and  only  If  there  Is  no 
scattering. 

Again,  considering  only  the  one-way  outgoing  wave,  we  deal  with  the 
solution 


♦  1k0  -  i*0/l  ♦  <n2  -  1)  ’  0. 


(10) 


Making  use  of  the  slowly  varying  property  of  n(r,z)  upon  r,  for  the  time 
being,  we  approximate  the  square-root  operator  by 


.••VO 
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(n2  -  1)  ♦  -X  i-y  si  ♦ 

k0  JZ 


(11) 


We  refer  to  this  as  the  small  angle  approximation.  Substituting  Eq.  (11) 
Into  Eq.  (10),  we  obtain  the  standard  PE,  that  Is, 


ur  -7ko(n2(r’z)  "  ')  u  +  ?Tuzz- 


The  derivation  Is  left  for  the  audience.  It  Is  Important  to  note  that  at 
this  point  the  standard  PE,  based  on  the  PE  approximation,  obeys  the 
following  limitations: 


LIMITATIONS 

1.  Farfield  Approximation,  k0r  ►  1. 

2.  n(i\z)  Slowly  Varying  in  r. 

3.  One-Way  Outgoing  Wave. 

4.  No  Scattering. 

5.  A  Particular  Square-Root  Approximation. 

VUGRAPH  4 
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Within  these  limitations,  the  standard  PE  Is  a  very  good  mathematical  model 
for  long  range,  low  frequency  propagation.  At  this  stage,  the  only  effective 
solution  algorithm  for  the  standard  PE  was  the  split-step  Fourier  algorithm 

3 

by  Tappert  and  Hardin. 

Before  I  mention  some  earlier  Important  developments,  I  want  to 
mention  the  'Workshop  on  Wave  Propagation  and  Underwater  Acoustics,*  held  In 
Mystic,  Connecticut,  In  November  1974.  This  workshop  was  not  limited  to  PE, 
but,  after  attending  this  workshop,  a  comprehensive  article  on  PE  was 

2 

written  by  Frederick  0.  Tappert  ("The  Parabolic  Approximation  Method,"  ), 
which  Is  a  chapter  of  the  book  cited  In  vugraph  5. 
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•  The  Parabolic  Equation  Approximation  in 

Wave  Propagation  &  Underwater  Acoustics 

F.  D.  Tappert 

ed.  J.  B.  Keller  and  J.  S.  Papadakis 

I  Lecture  Notes  in  Physics,  Springer-Verlag  (1977) 
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VUGRAPH  5 

To  this  date,  the  Tappert  article  Is  still  the  most  comprehensive  article  on 
PE. 
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In  1977,  there  was  another  workshop  In  Woods  Hole,  Massachusetts, 
where  Tappert  gave  another  paper  on  the  application  of  his  split-step 
algorithm.  Not  many  people  are  familiar  with  this  paper. 


Selected  Applications  of  the  Parabolic  Equation  Method 
in  Underwater  Acoustics 

F.  D.  Tappert 

International  Workshop  on  Low-Frequency  Propagation 
&  Noise,  Vol.  2,  Woods  Hole,  Mass.  pp.  155-194,(1974). 
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Since  1974,  the  Interest  In  the  PE  was  on  the  rise,  but  little  was 
done.  However,  there  were  a  number  of  notable  contributions. 


SrV* 
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The  first  practical  results  of  applying  the  PE  were  published  by  C.  W. 
Spofford.®  Then  a  few  papers  were  published  In  relation  to  the 
normal  mode  method  and  Its  application.  These  are 


if  f  t  f  fSS-' 


A  Synopsis  of  the  AESO  Workshop  on  Acoustic-Propagation 
Modeling  by  Non-Ray-Tracing  Techniques,  C.  W.  Spofford, 
A_E_SD_TN£3<)5_(1973j _ 

Eikonai  Approximation  and  the  Parabolic  Equation,  D.  R. 
Pajmer^^JVcoustr^oc  ._Am  ^f^2[,  Jl^-^£J_i976) _ 

Relation  Between  the  Solutions  of  the  Helmholtz  and 
Parabolic  Equation  for  Sound  Propagation,  J.  A.  DeSanto,  J. 
Acoj«tiSoa_Anv1 62i2)i_2Q5I297  j1977)_  _ _ 

A  Correction  to  the  Parabolic  Equation,  J.  A.  DeSanto,  J.  S. 
Perkins,  and  R.  B.  Baer,  J.  Acoust.  Soc.  Am.,  64(6),  1664-1666 
(1978J _ 

On  the  Parabolic  Approximation  to  the  Reduced  Wave 
Equation,  G.  A.  Kriegsmann  &  E.  W.  Larsen,  SIAM  J.  Appl. 
Math^  34^1)^201  ^04  ^1978)_ _ 

Helmholtz  Equation  as  an  Initial  Value  Problem  with 
Approximation  to  Acoustic  Propagation,  R.  M.  Fitzgerald, 

X  Acous t^_Soc^ Am^57(4] j/t 975[ _ 

Propagation  of  Normal  Mode  in  the  Parabolic 
Approximation,  S.  T.  McDaniel,  J.  Acoust.  Soc.  Am.,  57(2), 
30TS11_i1975) _ 

Parabolic  Approximation  for  Underwater  Sound  Propagation, 
S.  T.  McDaniel,  J.  Acoust.  Soc.  Am.,  58(6),  1178-1185  (1975) 
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Researchers  are  curious  and  Interested  In  the  relationship  between  the 
PE  and  the  Helmholtz  equation.  FI Itzgerald10  analyzed  the  PE  In  terms  of 
normal  mode  theory.  Dave  Palmer  6  observed  that  the  normal  mode  formulism 
was  difficult  In  ordering  the  geometric  optics  path-length  parameter  because 
of  mode-coupling.  The  PE  removes  this  difficulty. 


Then,  during  approximately  the  same  time  period,  a  few  Interesting 
developments  happened,  one  was  the  computer  code. 


The  Use  of  the  Parabolic  Equation  Method  in  Sound 
Propagation  Modeling, 

F.  B.  Jensen  and  H.  Krol 

SACLANTCEN  MEMO  SM-72  (1975) 

The  AESO  Parabolic  Equation  Model, 

H.  K.  Brock 

NORDA  TN-12  (1978) 

VUGRAPH  8 
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Jensen  had  a  PE  package  called  PAREQ  In  his  laboratory  performing 

14 

a  variety  of  research  and  applications;  so  too  did  Brock.  These  computer 
models  used  the  split-step  algorithm  to  solve  the  standard  PE.  G. 
Gartrell15  wrote  a  split-step  code  on  the  IBM  370/168. 

However,  In  solving  the  standard  PE,  a  number  of  users  found  a  phase 

O 

error.  DeSanto,  Perkins,  and  Baer  discussed  this  phase  error  and 
Introduced  a  correction  to  the  parabolic  approximation.  Notable  was  a 
technique  Introduced  by  Brock,  Buchal,  and  Spofford16  to  modify  the 
sound-speed  profile  to  Improve  the  accuracy  of  the  PE. 


Modifying  the  Sound-Speed  Profile  to  Improve  the  Accuracy  of 
the  Parabolic  Equation  Technique 

H.  K.  Brock,  R.  N.  Buchal,  and  C.  W.  Spofford 

J.  Acoust.  Soc.  Am.,  62(3),  pp.  543-552  (1977) 

VUGRAPH  9 

Interest  was  also  Increasing  In  the  application  of  the  PE  to  solve 
real  problems.  The  PE  models  In  various  laboratories  were  all  based  on  the 
use  of  the  split-step  Fourier  algorithm  with  an  artificial  bottom  treatment 


In  order  to  apply  the  PE  to  solve  real  problems,  the  model  would  be 
required  to  have  many  capabilities.  The  natural  question  Is  what  can  PE  do? 
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Can  PC  offer  these  capabilities.  These  questions  stirred  up  research  and 
development  Interest.  I  ask  the  question  In  a  different  way  —  what  can  we 
do  to  Improve  the  PE  capability? 

First,  to  Improve  the  PE  capability,  Lee,  Papadakls,  and 
Prleser1 ,f Initiated  the  numerical  solution  of  the  parabolic  wave 
equation  so  that  under  shallow  water  or  strong  bottom  Interaction 
environments  the  numerical  technique  can  handle  the  bottom  boundary 
condition. 


Numerical  Solution  of  the  Parabolic  Wave  Equation: 

An  Ordinary-Differential-Equation  Approach 

0.  Lee  and  J.  S.  Papadakis 

J.  Acoust.  Soc.  Am.,  68,  pp.  1482-1488  (1980) 

Generalized  Adams  Methods  for  Solving  Underwater  Wave 
Propagation  Problems 

D.  Lee  and  S.  Preiser 

J.  Comp.  &  Math,  with  Appls.,  7(2),  pp.  195-202  (1981) 
Finite-Difference  Solution  to  the  Parabolic  Wave  Equation 
D.  Lee,  G.  Botseas,  and  J.  S.  Papadakis 
J.  Acoust.  Soc.  Am.,  70(3),  pp.  795-800  (1981) 


VU6RAPH  10 

Lee  and  Papadakis  introduced  the  approach  of  using  an  ordinary 
differential  equation  (00E)  to  determine  a  special  kind  of  bottom  boundary 
(rigid)  within  the  framework  of  the  PE. 

This  bottom  boundary  treatment  was  Incorporated  Into  the  00E  and 

finite  difference  models.  At  that  time,  we  did  not  have  a  better 

Implementation  of  the  ODE  solution,  but  we  concentrated  our  efforts  on 

19 

developing  a  very  basic,  general  purpose  finite-difference  scheme,  which 

20 

could  be  Implemented  Into  computer  code.  This  scheme  Is  known  today  as 
the  Implicit  Finite  Difference  (IFO)  model  and  Is  an  Implicit  Crank-NIcolson 

9 


C«r 


method,  which  Is  unconditionally  stable.  The  IFO  model  is  used  quite  often 
to  solve  the  standard  PE. 

We  used  the  00E  model  to  solve  a  wedge  problem  consisting  of  range 
versus  propagation  loss.  The  results  turned  out  surprisingly  well. 


VUGRAPH  11 


After  we  showed  that  our  bottom  boundary  treatment  worked,  we  were  asked 
many  questions  pertaining  to  our  IFO  model's  capabilities.  One  of  the  first 
practical  questions  was  a  problem  proposed  by  H.  P.  Bucker.  He  asked:  Your 
IFO  model  can  handle  bottom  boundary  conditions  —  can  it  handle  the 
Interface  condition? 

The  problem  did  not  seem  to  be  a  particularly  difficult  one,  however, 
at  that  time,  no  existing  PE  code  could  do  It. 


DEPTH  (m) 


Handling  the  Interface  wasn't  difficult,  but  It  required  some  thought  as  to 

the  best  approach.  Thanks  to  many  valuable  discussions  with  Dr.  Suzanne  T. 

McDaniel,  we  worked  out  a  finite-difference  treatment  for  the  horizontal 
21 

Interface.  The  horizontal  Interface  development  Is  documented  In  the 
article  outlined  In  vugraph  14. 

A  Finite-Difference  Treatment  of  Interface  Conditions  for  the 
Parabolic  Wave  Equation:  The  Horizontal  Interface 

S.  T.  McDaniel  and  Ding  Lee 

J.  Acoust.  Soc.  Am.,  71(4),  pp.  855-858  (1982) 

VUGRAPH  14 

We  then  Incorporated  the  horizontal  Interface  conditions  Into  the  finite 
difference  code  and  ran  the  Bucker  problem  on  the  VAX  11/780  computer  at 
NUSC. 


VUGRAPH  15 


As  vugraph  15  shows,  without  the  Interface  treatment,  the  results  from  both 
IFO  and  split-step  do  not  agree  with  our  normal  mode  reference  solution 
results,  but  after  the  finite-difference  treatment,  the  results  are  In 

reasonable  agreement. 

12 


Naturally,  It  Is  logical  to  extend  the  finite-difference  technique  to 

22 

handle  the  Irregular  Interface  condition.  Dr.  McDaniel  and  1  have  worked 
this  out  and  I  shall  talk  about  this  a  little  later. 

While  the  capabilities  of  the  PE  were  being  developed,  use  of  the  PE 

model  was  Increasing.  Ourlng  this  period  of  development,  NOROA  sponsored  a 
workshop,  the  *N0ROA  Parabolic  Equation  Workshop. *23  Many  Interesting, 
realistic  problems  were  Introduced  at  the  workshop  and  a  comprehensive 
report  was  published. 


NORDA  Parabolic  Equation  Workshop 
James  A.  Oavla,  OeWayne  White,  and  Raymond  C.  Cavanagh 
NORDA  TN-143  (1981) 


VUGRAPH  16 

One  of  the  problems  Introduced  was  that  of  wide  angle  propagation. 


Frequency  =  290  Hz 
Water  depth  *  100  m 
Sound  apead  in  water  •  1900  nVa 
Oanaity  in  water  >  1.0  g/cm 
Oenaity  in  bottom  ■  1.2  g/cm 
Attenuation  in  water  «  0 
Bottom  attenuation  «  0.5  dB 
Bottom  sound  speed  a  1990  m/a 
Maximum  range  *  10  km 
Source  depth  ■  90.9  m 
Receiver  depth  «  99.5  m 
Number  ot  Modee  «  ii 


Sound  Speed  Profile 
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This  problem  was  solved  satisfactorily  by  both  the  normal  mode  and  the  fast 
field  program  (FFP).  We  then  used  the  FFP  solution  as  a  benchmark  reference 
solution. 
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VU6RAPH  18 

Note,  In  vugraph  18,  that  the  dashed-dotted  line  was  the  solution  produced 
by  the  standard  PE  (both  the  split-step  and  the  IFO),  the  disagreement  Is 
clear.  It  was  due  to  the  size  of  the  angle  of  propagation.  An  Important 
capability  required  to  produce  agreement  Is  the  wide  angle  capability. 


Let  me  briefly  describe  the  mathematical  development  of  the  wide  angle 
capability.  Recall  the  one-way  outgoing  wave  equation,  that  Is, 


(n2  -  1)  ♦  -j2-?  )u  *  °. 
ko  iz  / 


where  we  made  the  approximation 


Jl  *  [ n2  -  1)  ♦  -l  *“?  » 1  ♦  i  (n2  -  1)  ♦ 


to  obtain  the  standard  PE 


ur  -7ko(n2(r-z)  -  l)  u  uzz' 


One-way  outgoing  wave:  iK°  "  ik°  y1  +  ~ 

J 1  +  W  -  1>  +  1  =  1  +  1  f(n»  -  1)  +  1_ 

’  '  **,3*  7  [  •£ 


=  0 


IL 

!32! 


ur  =-|K0<nJ(r,z)  -  b u  + 
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To  be  general,  we  approximate  the  square-root  operator  by  a  rational 
function  approximation  as  follows: 


*o 


If  p  -  1/2  and  q  -  0,  we  see  clearly  It  reduces  to  the  small  angle 
approximation,  thus,  resulting  In  the  standard  PC. 


..[(A,.,) 


1  ♦  qfc^r.z)  -  l)+  ^  f 
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(12) 


For  a  special  selection,  p  »  3/4  and  q  -  1/4,  the  right-hand  side  of 
Eq.  (12)  becomes 


VlT^r.z)  -  l)*^ 


—7  — T* 


1*1 

(n2(r,z)  - 

')*44l 

• 

k.  az 

o  . 

>♦) 

(n2(r,z)  - 

l)  ~7  ~~7 

’  *o  ”  J 

• 

(13) 


which  Is  the  rational  function  approximation  of  the  square-root  operator  by 
25 

John  F.  Claerbout.  We  choose  to  keep  p  and  q  arbitrary  so  that  we  can 
determine  p  and  q  to  suit  our  needs. 


If  we  use  Cq.  (12)  for  the  square-root  operator  and  substitute  It  Into 
the  one-way  outgoing  wave  equation,  we  find  a  pseudopartlal  differential 
equation,  that  Is, 


! 


(14) 


which  Is  the  Wide  Angle  Wave  Equation. 
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A  Rational  Function  Approximation 
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Jon  F.  Claerbout  approximation 


Wida  Angle  Wave  Equation 
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We  tend  to  refer  to  It  as  the  Wide  Angle  PE,  but  In  a  sense  this  Is 
deceptive.  The  PE  Is  our  key  equation.  When  p  -  1/2  and  q  *  0,  the  standard 
PE  Is  a  special  case.  Therefore,  when  we  refer  to  the  Wide  Angle  PE,  we  are 
actually  referring  to  the  pseudopartlal  differential  equation  (14). 

With  the  wide  angle  capability,  you  see  that  the  IFD  produced 
agreeable  results  when  compared  with  the  benchmark  FFP  solution. 

- Fast  field  (FFP  -Exact") 


- Standard  PE 

....  Implicit  finite  difference- 
(Claerbout  equation! 


5.0  5.5  6.0  6.5  7.0  7.5  8.0  8.5  9.0  9.5  10.0 

RANGE  (kmf 

Wide  Angle  Solution  Comparison 
VUGRAPH  21 
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There  have  been  a  number  of  authors  who  have  contributed  to  the 
theoretical  development  of  the  wide  angle  capability. 


Fundamentals  of  Geophysical  Data  Processing  with 
Applications  to  Petroleum  Prospecting 

Jon  F.  Claerbout 
McGraw-Hill  (1976) 


High  Angle  PE 
Robert  R.  Greene 

NORDA  Parabolic  Equation  Workshop  (1981) 

IFD:  Wide  Angle  Capability 

George  Botseas,  Olng  Lee,  and  Kenneth  E.  Gilbert 

NUSC  TR#8905  (1983) 

Extension  of  the  Parabolic  Equation  Model  for  High-Angle 
Bottom-Interacting  Paths 

L.  B.  Dozier  and  C.  W.  Spofford 

SAI  Technical  Report  SAI-78-712-WA  (1977) 
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Greene  used  the  term  "rational  parabolic"  as  opposed  to  "parabolic."  He  also 

applied  the  rational  function  approximation  to  the  square-root  operator. 

This  rational  function  approximate  technique  was  applied  earlier  by 
25 

Claerbout. 

A  number  of  other  people  also  examined  and  developed  the  wide  angle 

capability.  At  NUSC,  we  Incorporated  this  Into  our  IFO  code  and  a 

24 

comprehensive  report  was  published.  In  the  seismology  field,  Berkhout 

28 

used  the  continued  fraction  to  approximate  the  square  root  operator. 


Wave  Field  Extrapolation  Technique  in  Seismic 
Migration 

A.  J.  Berkhout 

Geophysics,  Vol.  46,  #12,  pp.  1838-1656  (1981) 
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It  Is  Interesting  to  note  that  Berkhout's  first  order  continued  fraction 
coincides  exactly  with  Tappert's  standard  small  angle  PE.  His  second  order 
continued  fraction  coincides  with  the  Claerbout  approximation.  Roughly 
speaking,  the  small  angle  PE  can  accommodate  propagation  angles  up  to  15*. 
and  wide  angle  PE  (IFO)  can  accommodate  propagation  angles  up  to  40*. 
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David  Thomson  was  the  first  to  apply  the  Split-step  algorithm  to 
handle  the  wide  angle. 


A  Wide  Angle  Split-Step  Algorithm  for  the  Parabolic 
Equation 

David  J.  Thomson  and  N.  R.  Chapman 
J.  Acoust.  Soc.  Am.,  74(6),  pp.  1848-1854  (1983) 
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Recently,  one  of  my  colleagues,  Donald  St.  Mary  of  the  University  of 
Massachusetts,  developed  a  very  high  angle  PE  using  a  higher  order  rational 
approximation  that  can  accommodate  a  propagation  angle  greater  than  45*. 
This  Is  good  news  In  the  area  of  shear  wave  propagation. 


Formulation  and  Discretization  of  a  Very  Wide  Angle 
Parabolic  Equation 

Donald  F.  St.  Mary  and  Ding  Lee  (1984) 
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All  the  authors  I  have  mentioned  thus  far  have  made  contributions  to 

the  wide  angle  capability  In  this  period  of  time.  Two  researchers,  Estes  and 
31 

Fain,  started  earlier  examining  the  wide  angle  formulation  using  the 
Taylor  series  expansion. 
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Numerical  Technique  for  Computing  the  Wide  Angle 
Acoustic  Field  In  an  Ocean  with  Range-Dependant 
Velocity  Profiles 

L  E.  Estes  and  G.  Fain 

J.  Acoust.  Soc.  Am.,  62(1),  pp.  38-43  (1977) 


VUGRAPH  26 


Continuous  use  of  the  PE  for  research  and  application  purposes  showed 

reasonable  success,  however,  under  unusual  environments.  It  was  not  fully 

developed  to  handle  everything.  SACLANT  and  NOROA  called  our  attention  to 

the  fact  that  In  the  Irregular  sloping  Interface  situation.  If  the  depth 

partition  points  do  not  fall  on  the  Interface  boundary.  Inaccuracy  will 

occur.  In  numerical  analysis  language,  this  Is  called  numerical  reflections. 

The  present  computer  code  cannot  handle  the  situation  without  some 

modification  because  the  Irregular  Interface  conditions  treatment  Is  not 

Included  In  the  present  code.  In  addition  to  the  progress  made  by  McDaniel 
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and  me,  recent  progress  has  been  made  by  Jaeger  to  treat  the  Interface 
condition  by  using  the  Irregular  Interface  condition  developed  by  McDaniel 
and  me. 


A  Computer  Program  for  Solving  the  Parabolic  Equation 
Using  an  Implicit  Finite-Difference  Solution  Method 
Incorporating  Exact  Interface  Conditions 

Larry  Ernest  Jaeger 

Naval  Postgraduate  School,  MS  Thesis  (1983) 

A  Finite-Difference  Treatment  of  Interface  Conditions  for  the 
Parabolic  Wave  Equation:  The  Irregular  Interface 

Ding  Lee  and  S.  T.  McDaniel 

J.  Acoust.  Soc.  Am.,  73(5),  pp.  1441-1447  (1983) 
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Independently  from  the  direct  application  of  the  irregular  Interface 

33 

condition,  Jules  deG  Grlbble  made  an  Important  Improvement  on  the  finite 
difference  model. 

Jules  allows  variable  mesh  spacing  and  deals  with  Interfaces  that  do 
not  lie  on  a  mesh  point.  He  extended  the  finite  difference  treatment  of 
Interfaces  to  the  parabolic  wave  equation.  Jules1  main  effort  was  centered 
at  the  employment  of  variable  vertical  mesh  size  In  conjunction  with  a  high 
level  system  of  00£  solver. 


Extending  the  Finite  Difference  Treatment  of  Interfaces 
When  Using  the  Parabolic  Wave  Equation 

Jules  deGribble 

J.  Acoust.  Soc.  Am.,  76(1),  pp.  217-221  (1984) 
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Up  to  this  point.  It  seems  that  there  Is  enough  PE  capability  for  research 
studies  as  well  as  for  applications.  Three  types  of  PE  models  exist. 


EXISTING  PARABOLIC  EQUATION  MODELS 

1.  Split-Step  Fourier  Algorithm  Model 

2.  Implicit  Finite-Difference  Model 

3.  Ordinary-Differential-Equation  Model 
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Little  has  been  done  to  compare  the  available  computer  models.  The 

34 

only  published  literature  was  given  by  Kewley  who  discussed  practical 
solutions  of  the  PE  model  for  underwater  acoustic  wave  propagation. 


20 


v- 


Practical  Solutions  of  the  Parabolic  Equation  Model  for 
Underwater  Acoustic  Wave  Propagation 
in  Computational  Techniques  &  Applications  Conference 

0.  J.  Kewley,  L  T.  Sin  Fai  Lamn,  and  G.  Gartreli 
ed.  J.  Noye,  Sydrey 

North-Hoi  land,  Amsterdam  (1983) 
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Among  these  three  models,  IFO  Is  the  more  general  purpose.  The  ODE 
solution  has  great  potential,  but  Is  not  yet  fully  developed.  Applications 
of  the  PE  model  Indicates  that  It  Is  doing  well  as  a  research  code, 
especially  the  IFD  code,  because  of  Its  accurate  computation  of  the  wave 
field.  However,  a  number  of  users  are  actually  using  these  codes  for 
applications. 

In  solving  either  the  small  angle  PE  or  the  wide  angle  wave  equation, 

we  deal  with  one  Input  parameter,  the  reference  wavenumber  kQ.  Application 

of  the  above  equations  to  solve  real  problems  requires  a  clever  selection  of 

the  reference  wavenumber  ko>  This  selection  was  Ignored  by  the  model 

developers  and  the  users.  For  the  user.  In  practice,  he  does  not  usually 

have  the  knowledge  to  select  the  best  k  ;  therefore,  users  have,  in  most 

0 

cases.  Ignored  the  selection  of  kQ.  The  model  developer  built  In  the 
automatic  selection  of  kQ  according  to  physical  experiences.  Inappropriate 
kQ  will  lead  to  an  evident  phase  error.  Pierce  recently  reemphasized  the 
Importance  of  kQ  selection  and  Introduced  a  formula  to  determine  the  range 
of  kg  based  on  the  Rayleigh  quotient.  Some  numerical  experiments  have  been 
carried  out  at  the  Naval  Underwater  Systems  Center,  New  London  Laboratory; 
results  show  some  phase  shift  effects  dependent  on  kQ  variations.  Tappert 
and  Lee  joined  Pierce  In  studying  the  natural  selection  of  k0- 
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The  Natural  Reference  Wavenumber  for  Parabolic 
Approximations  in  Ocean  Acoustics 

Allan  D.  Pierce 

J.  Comp.  &  Math,  with  Appls.  (To  be  published)  (1984) 
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It  is  interesting  to  note  that  L.  Ngheim-Phu  and  Tappert  developed 
a  PE  box  that  can  be  put  on  shipboard  for  fast  prediction. 


A  High-speed.  Compact,  and  Interactive  Parabolic  Equation 
Solution  GENerator  System  (PESOGEN) 


Lan  Ngheim-Phu  and  F.  D.  Tappert 
J.  Acoust.  Soc.  Am.,  75(S1),  p.s26  (1984) 
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The  interest  in  the  PE  solution  has  continuously  increased.  An 
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interesting  development  was  the  Yale  University  Workshop.  Three  days 
were  devoted  to  discussions  by  Invited  speakers  about  the  state-of-the-art 
of  computational  ocean  acoustics. 


COMPUTATIONAL  OCEAN  ACOUSTICS  WORKSHOP 

YALE  UNIVERSITY 

Sponsors:  ONR  (Math.  Dept.) 

Yale  U.  (Computer  Science  Dept.) 

NUSC  (Independent  Research) 

Special  Issue:  J.  Comp.  &  Math,  with  Applications 


Pergamon  Press 

Editors:  Martin  H.  Schultz  and  Ding  Lee 
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Range  Refraction  Corrected  PE.  (ref  40) 

Density  Effects  in  PE.(ref  39) 

Time  Domain  PE. 

High  Frequency  PE  (hybrid  ray-PE).  (ref  39) 
Backscattering  (in  range)  PE.(ref  41) 

Rough  Boundaries  in  PE.(ref  42) 

3-Dimensional  PDE.(ref  43) 

Stable  Explicit  Schemes  for  Solving  3-Dimensional 
(2-Dimensional)  PE 
Shear  Waves  in  PE. 

Match  Interface  Conditions  Between  Fluid  and  Elastic 
Media. 

Large  Angle  PE. 

Currents  in  PE. 

Ooppler  Effects  in  PE. 

Moving  Sources  and  Receivers,  or  Time-Variables. 

The  Application  of  Multi-Array  Processors  to  Solve  the 
Parabolic  Wave  Equation 

Reciprocity  in  the  Time  Domain  and  the  PE  Method. 
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As  for  the  computation  of  long  range  propagation.  It  Is  on  a  very 
large  scale.  Ideally,  we  would  like  to  have  In  principle 


LARGE  SCALE  COMPUTATIONS 


1.  Accuracy 


2.  Speed 


3.  Minimal  computer  storage  requirements 
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Let  me  use  the  IFO  as  a  simple  example.  Using  the  IFD  to  solve  the  parabolic 
wave  equation,  either  small  angle  or  wide  angle,  requires  the  solution  of  a 
system  of  equations  of  the  form 


Aun  + 1  =  Bun  +  un  +  un  + 1 
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Since  both  matrices  A  and  B  are  trldlagonal,  a  special  recursive  formula  Is 
applied  to  solve  the  system  at  minimal  cost.  Not  counting  the  overhead,  the 
trldlagonal  solver  requires  6N  operations.  The  solution  looks  so  simple  and 
economical  that  It  does  not  seem  to  require  Intensive  computations  even  for 
long  range  propagation  problems.  However,  this  Is  not  generally  true.  Let  us 
recall  our  earlier  solution  to  the  wedge  problem  of  a  rigid  bottom  boundary 
condition.  Our  early  solution  to  this  problem  was  applying  the  Generalized 
Adams  00E  method,  which  produced  the  accurate  solution. 
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For  the  sloping  bottom  boundary,  In  order  to  solve  the  system 
satisfactorily,  we  adopted  a  variable  dimension  procedure.  We  started  at 
approximately  348  m  for  a  0.15  m  depth  Increment  and  ended  up  to  solve  an 
Initial  system  of  equations  of  order  400.  Due  to  the  memory  storage  limit, 
calculations  were  not  allowed  to  go  up  to  a  maximum  range  of  10  km,  which 
means  we  have  to  solve  a  system  of  equations  of  order  5740.  Recently,  we 
Improved  the  technique  and  relaxed  the  storage  requirement,  thus.  Increased 
the  speed. 
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The  parabolic  wave  equation  can  be  solved  by  a  first  order  Generalized 
44 


Adams  method. 


PE 

ur  =  a(k0,r,z)u  +  tyko.r.z)^ 

ODE 

Ur  =  A(k0,r,z)u  +  g(k0,r,z,u0) 

GAB 

um  + 1  =  eAhun  +  h(Ah)  -  1(eAh  -  |)gn 

EXP 

8Ah=(|_Ah)-1 

Improved 

(1  -  Ah)un  + 1  =  un  +  hgn 

GAB 
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Since  Matrix  A  Increases  dimension  at  every  range  step,  e  has  to  be 

3 

updated  at  every  range  Increment.  It  costs  N  operations  for  the  Inversion 
2 

and  cost  N  memory  storages  because  we  calculate  the  matrix  exponential  by 
the  rational  Pade  approximation,  which  requires  the  Inversion  of  an  NxN 
matrix.  Even  though  the  Matrix  A  Is  tcldlagonal,  the  Inverse  destroys  the 
trldlagonal  property  and  fills  up  the  storage.  You  can  see  the  cost.  By  a 
careful  study,  for  a  reasonable  choice  of  the  step  size,  we  can  calculate 
the  new  wave  field  by  avoiding  this  expensive  matrix  exponential  calculation 
and  solving  a  system  of  equations. 

• 

3 

This  system  Is  again  trldlagonal,  we  reduce  N  operations  to  6N 

2 

operations;  moreover,  we  reduce  a  N  memory  storage  requirement  to  only  3N 
storages.  By  coincidence,  this  problem  can  also  be  solved  nicely  by  the  IFD. 
Now,  one  can  again,  by  taking  advantage  of  using  hardware,  Improve  the 
efficiency  of  the  computation.  You  shall  hear  this  from  other  talks. 


■v.vv.vv 
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CONCLUSIONS 

Now  you  have  heard  of  the  many  Important  contributions  toward  the 
application  of  PE  approximation  In  solving  ocean  acoustic  problems  and  have 
also  heard  of  some  Interesting  PE  developments.  You  may  note  that  there  are 
many  not  yet  well -developed  PE  approximations,  many  more  than  existing  PE 
capabilities.  However,  I  would  like  to  call  your  attention  to  the  fact  that 
NOT  every  problem  can  be  solved  by  the  PE  approximation  In  an  efficient 
manner.  If  the  physical  conditions  fall  within  the  limitation  of  the  PE 
approximation,  PE  approximation  Is  an  efficient  method  to  apply. 


NOT  every  problem  can  be  solved  by 
the  Parabolic  Equation  Approximation 

VU6RAPH  41 

I  welcome  your  comments  and  participation  In  solving  these  problems  and 
applying  your  solutions  to  realistic  problems. 


Thank  you. 
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